Shear viscosity of a superfluid Fermi gas in the unitarity limit 
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We compute the shear viscosity of a superfluid atomic Fermi gas in the unitarity limit. The unitarity limit is 
characterized by a divergent scattering length between the atoms, and it has been argued that this will result in 
a very small viscosity. We show that in the low temperature T limit the shear viscosity scales as /T^ , where 
the universal parameter ^ relates the chemical potential and the Fermi energy, fj = t,Zf. Combined with the high 
temperature expansions of the viscosity our results suggest that the viscosity has a minimum near the critical 
temperature Tc. A naive extrapolation indicates that the minimum value of the ratio of viscosity over entropy 
density is within a factor of ~ 5 of the proposed bound T|/.s > h/i^nkB). 



I. INTRODUCTION 

Shear viscosity r| can be defined as the shearing force F 
per unit area A per unit velocity gradient in a laminar flow. 
For a flow in x-direction, with a velocity gradient VyVx in the 
y-direction 
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Viscosity relates the rate of momentum transfer to the veloc- 
ity gradient. For dilute gases the microscopic mechanism for 
momentum transfer is provided by atomic collisions. This 
mechanism becomes more efficient as the mean free path gets 
larger because in that case the atoms travel larger distances be- 
tween collisions and transfer momenta between laminar layers 
of more disparate flow velocities. Thus viscosity r| is expected 
to be inversely proportional to the collision cross section o. 
This leads to the question of whether there is a fundamental 
limit to how small the viscosity can get as the strength of the 
interaction is increased. Stated differently, we would like to 
determine the shear viscosity of the most "perfect" fluid. 

There is an old argument that suggests that quantum me- 
chanics places a lower limit on the shear viscosity Ul]. A 
rough estimate of the viscosity is provided by r| np'k, where 
n is the number density, p is the average momentum, and A, the 
mean free-path. Heisenberg's uncertainty principle requires 
pX> h and the kinematic viscosity r[/n > h. For relativis- 
tic systems particle number is not conserved and it is more 
natural to consider r|/i, where s is the entropy density. As 
long as the entropy per particle is of the order kg we expect 
r|/i > h/kB. 

A new perspective on this idea is provided by a calculation, 
based on the AdS/CFT correspondence, of r|/i in the strong 
coupling limit of = 4 super-symmetric Yang Mills theory 
||2tl. This calculation gives r\/s^h/ (471^^), a value that is also 
obtained in other strongly coupled field theories that have a 
gravity dual. It is also known that the leading order correction 
to the limit of infinite coupling increases r|/i. This has led 
to the conjecture that the strong coupling result is a universal 



lower bound for all fluids ^ : 
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Liquid Helium comes to within an order of magnitude of the 
bound, and values r[/s ~ (0.1 — 0.5)h/kB have been reported 
for the quark gluon plasma produced at RHIC 10, Ht]. There 
are suggestions in the literature that counter examples can be 
found by considering non-relativistic systems for which the 
entropy per particle is very large but currently no fluid 

that violates the bound is experimentally known. 

An interesting system to study in this context is a cold 
atomic gas near a Feshbach resonance |^ [T^ [TH [T2I1 . In 
^Li and '^''K gases, there exist hyperfine channels that sup- 
port bound states. The magnetic moment of the bound state 
in these channels is different from the sum of the magnetic 
moments of the atoms that make the bound state. This al- 
lows one to use an external magnetic field to move the bound 
state energy relative to the continuum states, effectively mak- 
ing the bound state arbitrarily shallow. In terms of scattering 
theory, a shallow bound state corresponds to a large scattering 
length. At the Feshbach resonance, the atomic cross section 
is only limited by unitarity a{k) ~ l/k^. The unitarity gas 
interaction is characterized by a divergent two-body scatter- 
ing length |fl| ^ 00 and a natural sized range r ^ lA. Even 
for a dilute gas with density n ^ r^^, the unitarity gas with 
|fl| ^ 00 is a strongly interacting system. In fact it is the most 
strongly interacting non-relativistic system known, with a di- 
verging two-body collision cross section a{k = 0) ~ a ±0°. 

The aim in this work is to improve the understanding of 
transport properties of the cold unitarity gas by performing a 
systematic calculation of the shear viscosity in the low tem- 
perature superfluid phase. Combined with known results in 
the high temperature limit ifisll these results provide an es- 
timate of the minimum viscosity. In the superfluid phase 
Cooper pairs break the U{1) symmetry associated with the 
conservation of particle number. This implies that there is a 
Nambu-Goldstone boson, the phonon. At temperatures T be- 
low the critical temperature Tc for superfluidity, phonons dom- 
inate thermodynamic and transport properties of the system. 

The paper is organized as follows. In Section|II]we present 
the basic equations relating the shear viscosity to the phonon 
collision operator. The phonon interaction is derived in Sec- 
tion |III1 followed by a variational calculation of the viscosity 
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in Section |IV] A discussion of the result is presented. The 
discussion closely parallels the calculation of the viscosity in 
hquid Helium 01 H ' 
dense quark matter 1116 
tion[Vl 



and, in particular, the CFL phase of 
We end with the conclusions in Sec- 



n. TRANSPORT EQUATION AND VISCOSITY 

Viscosity as defined in Eq. ([T|i is related to internal stresses 
in a fluid. A more convenient definition is provided by the 
stress-energy tensor 7) ,- of an almost ideal fluid. Close to equi- 
librium it can be expanded in derivatives of the flow velocity 



T,j=iP + e)ViVj-P5ij + 5Tij , 

bT,j = -^ {ViVj + ViVj - ^bijV ■v) + - 



(3) 



where we only kept the traceless part of 57;^. The trace of 
dTij is related to bulk viscosity. The ideal fluid part of Tjj is 
related to the thermodynamic variables pressure P and energy 
density e. In the superfluid phase the long distance fluctua- 
tions of the order parameters and of the conserved quantities 
are described by the two-fluid hydrodynamics. The two com- 
ponents are a non-viscous superfluid, and a viscous normal 
fluid. The stress-energy tensor of the normal fluid is given by 
Eq. (O, where V,- is now the velocity of the normal fluid. 

If the normal fluid is composed of weakly interacting quasi- 
particles the stress-energy tensor and the viscosity can be 
computed using kinetic theory. In the unitarity Fermi gas at 
very low temperature the quasi-particles are the phonons. The 
stress-energy tensor is given by IitIi 



T.J = 



d^P PiPj 



fp, 



(4) 



where fp is the distribution function of the phonons with 
speed V, momenta pi and energy Ep. Close to the equilibrium 

fp — fi'^^ + ^/p' where /p"^ is the Bose-Einstein distribution 
and 8/p is a small departure from equilibrium. Small fluctua- 
tions can be parameterized in terms of departures of the ther- 
modynamics variables r,/j,y,- from equilibrium, e.g. 5/p ^ 

Tdrfi^^ - +fp°^)/T. This motivates the definition 



8/p 



-X{p)fi7^\l +fp'^)/T in terms of the unknown func- 



tion x{p). To project onto the shear stress, one uses the ansatz 

Xip) ^gip)iP.Pj^l5ijP')(y.Vj + '7,iV,-^5,jV-V), (5) 

where only the traceless projection on the momenta p is rele- 
vant. Thus close to the equilibrium one can write 



d^P PiPj 

(27t)3 Ep 
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x(v,y, + v,y,--5,;V-v). 



This determines the shear viscosity in terms of the function 

8{P), 
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The equation of motion for g{p) is derived using the Boltz- 
mann equation 



dfp ^dfp 
dt dt 



'V-Vfp + F-Vpfp^C[fp], 



(8) 



relating the rate of change of the distribution function fp to 
the collision operator C[fp]. In the absence of external force 
we take F = 0. The left hand side of the relation Eq. (O can 
be simplified further iITtIi to write; 



"Jp _ Jp , AQ)u ^ s 2\ 

x{V,Vj + VjVi-^8ijV-V), 



(9) 



where only the contribution relevant for shear viscosity was 
retained in the Linear Response Approximation, leading order 
in the small deviation from equilibrium. 

Two types of contributions to the collision term C[fp] are 
typically considered: (a) binary 2^2 collisions in which the 
number of particles is conserved, and (b) 1 — > 2 "splitting" 
processes in which the number of particles is not conserved. 
These processes are shown in Fig. [T] We will show in Sec- 
tion |III]that splitting processes do not contribute to shear vis- 
cosity at leading order in the low temperature approximation. 

The 2^2 collision integral is given by 



1 r d^k 
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d^'k 



d'p' 



lf2Ek {2K)^2Ek, {2%f2Ep, 

X (27t)45W(p + yt-p'-/t')|M|2Z)2«2, 
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where Di^i contains the distribution functions and |M| is the 
2^2 scattering amplitude shown in Fig.[T] The distribution 
functions are linearized in the small deviations from the equi- 



librium distribution, D2^2 ~ D 



(0) 



bD2^2- We find: 
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where we have used the equilibrium relation 

/f/.i^i +/r)(i +4?') =(1 +/r)(i +/r 



This relation ensures that D 
equilibrium. 
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FIG. 1 : The first two diagrams contribute to the 2 2 process, and 
the last four diagrams contribute to the 1^2 processes. Only the 
leading order contribution to the shear viscosity T| from the vertices 
are included. 



Using the ansatz in Eq. (|5]l for x{p), we get 

C2^2[fp 



^+fl^^ [y , /I , ^(0)^^(0) .(0) 

2^^j J ^k-k'p'{^+Jk )Jk' Jp' 



(13) 



X [g{p)p.j+8ik)k,j~8ik')k'u^8ip')Pij] Vij 
=Fij[g{p)]Vij, 

where we have defined the linearized collision operator 
Fij[g{p)]. We have also defined 



3 „/ 
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{InyiEk {InyiEk, {InflEp, 
x{2%)W\p + k-k' -p')\'M\^, 



(14) 
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exp{Ep/T) - 1 



Using Eqs. (O, (|9]l and ( fT3] l. the linearized Boltzmann equa- 
tion can be written as 



FijHp)] =^ 



2pT 



Pij- 



(15) 



This result can be used to rewrite the relation for the viscosity 
in Eq. (|7| as 



^=5 



d^p 



Pij8ip)Fij[g{p)h 



(16) 



which will be useful later. We used the linear dispersion rela- 
tion Ep = vp above, sufficient for the calculation as shown in 
the next section. To complete the calculation of the solution 
to the collision equation we need to calculate the scattering 
amplitude ~ |M| which we will turn to now. 



III. PHONON CROSS SECTION 

The phonon interaction for the unitarity gas in the super- 
fluid phase can be derived from Galilean and gauge invari- 
ance iflsllTgll . Consider a microscopic Lagrangian for the uni- 
tarity Fermi gas 



2m 



-P 



(17) 



where \|/ is two component spinor, m is the mass of the 
Fermion, 02 is the anti-symmetric Pauli matrix, and Co is an 
interaction strength that can be tuned to achieve infinite scat- 
tering length. This Lagrangian is invariant under Galilean 
transformations, and under the gauge transformation \\t 
e"?W\|/ where the fictitious gauge field Ay Av — d^q is de- 
fined as Av = {fJ,0). We work in units where h = I ^ c = ks- 
We require that the effective theory for the phonon field (|) 
shares the symmetries of the microscopic Lagrangian. This 
implies that the effective Lagrangian L,^ is a function of 

X = A'-5oC^-(V(^)V(2'«), (18) 

and its derivatives lfl9ll20ll . The functional dependence on % 
is further restricted by the observation that the effective action 
at its minimum r(x ^fS) = T j d^xLi^, for constant classical 
field 8v(|) = 0, is equal to the pressure of the unitarity gas. 
In the limit |al — > 00, r = which is nearly realized in cold 
atomic traps lH H, [ii [H d, the unitarity gas is a scale in- 
variant system. This implies that, up to a numerical constant, 
the pressure P has to be equal to that of the free system. We 
write 



,5/2 



(19) 



where the universal constant ^ is sometimes called the Bertsch 
parameter in the nuclear physics community. Eq. ( fT9l ) implies 
^ = , where e/r = kp/(2rn), kf = (SJt^n)'/-', and n is the 
number density. We conclude that 
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2n 5/2 



2m 



where . . . corresponds to terms with derivatives of %. We 
can bring the kinetic term into the canonical form via a field 
rescaling (|) 7i^^''^(|)/[(m^/j)^'''*2'/'*]. We find expanding in 
derivatives of the phonon field, ignoring total derivatives of 
the dynamical field (|) and constants independent of (|), 



(21) 



a 



-2«'[(3oC^) 



9v28o^(V0)2] 

* + 18v2(8o^)2(Vc^)2 - 27v\V(^)*] 



where a = 7tv^/2^3/4^j-3i/4g^2-j ^j^^ Nambu-Goldstone 
boson speed is = 2^/ {3m). 
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FIG. 3: Leading order contributions to tlie phonon self-energy cor- 
rection. The tadpole does not contribute an imaginary part. 



FIG. 2: Leading order contributions to the binary collisions. 



The determination of ^ is a non-perturbative many-body 
problem, and there are no exact analytical calculation avail- 
able. Numerical calculations using fixed node Green's Func- 
tion Monte Carlo ll2lll22ll23ll or Euclidean lattice calculations 
iMmEilll find ^ - 0.3 - 0.4. Our final result depends on 
this single universal number t,. 

We can estimate the sizes of the different terms in the La- 
grangian as follows: for the kinetic term to contribute to the 
generating functional its contribution should be 0(1) other- 
wise it will be damped in the exponential. Time deriva- 
tives scale as 8o ^ T, spatial derivatives as 8,- ~ T/v, and 
the volume integral scales as d'^x ^ /T^. This implies that 
(|) ~ r/v^/^. We observe that the magnitude of the phonon 
self coupling relative to the kinetic term scales as a(3o(|)) ~ 
^l'^{T /nY, a small correction for T <^ /u. Note that for a 
strongly interacting unitarity gas Tc ~ (0.29±0.02)7> w 0.7/j 
ll28ll for ^ = 0.4, which implies that Tc is of the order /j. 

The Lagrangian in Eq. Uli describes the leading order 
phonon interaction for the processes shown in Fig. [T] At this 
order, the phonon dispersion relation is linear with Ep = v\p\. 
Consequently, the splitting processes 1 ^ 2 are collinear and 
cannot contribute to the shear viscosity. 



Binary Collisions 

The leading order contribution to the binary collision pro- 
cesses in Fig. [T] are shown in Fig. |2] The contribution 
of the four-phonon contact term to the scattering amplitude 



32^(4 



kpk'{3 [cosy(fc — 6fccos9 — 3/:') 



+k + p+ [p -3k) cose' + cose(yt + /9 - 6/9Cose')] 

+ 3 [3 cos e - cos e' + cos y(6 cos e' - 1 ) - 1 ] /t' 

+ 3(cosy+cose + cose')p' + p'}, (22) 

where we have used p + k = p' + k' and defined p-k = cos 0, 
p ■ k' = cos 9' and k ■ k' — cosy. We also assumed that the 
phonons are on-shell and that the dispersion relation is linear, 
= ^'\p\- Factors of 1/2 from Bose symmetry have been 
included in the amplitudes. 



If the phonon dispersion relation is linear the s, t and u- 
channel phonon exchange amplitudes diverge in the collinear 
limit. This corresponds to sub-sequent collinear splitting and 
joining processes with an on-shell propagator in between. The 
collinear processes should not contribute to the shear viscos- 
ity, but the numerical evaluation of collision integrals is more 
stable if the infrared divergence due to the on-shell propagator 
is regularized by including the thermal damping of the phonon 
propagator For this purpose we compute the imaginary part 
of the self-energy correction T.{p) to the phonon propagator. 

There are two self-energy diagrams at O(a^), Fig. [3] The 
tadpole graph does not generate an imaginary part and we only 
compute the first diagram. We find 



yf ^ tf^T V 



d^k 



1 



(23) 



16x/3a(4 „£^J (27t)3(o2+£2 

[po{2P ■K-K^)+ h){P^ - IP ■ K)] ^ 
{~ipo + Q+-(£>„)^+Ek-p 

The four- vector products are defined as P K= pok n — 9v^p ■ k, 
P^ ~ Pq — 9v^p^. This can be computed following 111611 and we 
find 



^PO,P) 



^2^,3^3/2 ^ ^ ^3^ 

16\/3/ 

po + ;0+ - siEi, - S2Ep_k 
[po{2P ■K-K^)+ ko{P^ - 2P ■ K)] 



(24) 



The imaginary part of 'L{pQ,p) arises from the pole terms in 
the propagator. Analytic expressions for lm'L{po,p) can be 
found in Appendix|A] For very time-like \po\ ^ \p\ external 
momenta 



Imi:(po,p) 



3V37tj,3/2 6 

exp(^) + l 
"exp(^)-l 



exp(# 
exp(§) 



1 



W—T®iPO 

TT I ^ 



(25) 



@hpo) 



and for space-like \pQ \ < \'\p\ external momenta with v\p\<^T 

lmJ:ipo,p) « ^^3/2j,4/i0(^2 2_ 2) (26) 

5^*v p 
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For the calculation these limiting forms provide sufficiently 
accurate representations of the exact one-loop expression in 
Eq. (l24l) . We define the dressed phonon propagator 



iG{po,p) 



Pq-v^P^ + ilmY.{pQ,p) 



(27) 



We can now collect the regularized s, t and M-channel phonon 
exchange amplitudes. The i-channel amplitude is 

iM, = -/ ^^^^ (p + kfG{po + ko,p + k) (28) 

X [4-v^pk-P-K] [Av^p'k' -P' -K'] . 

The t and M-channel amplitudes follow from crossing symme- 
try. /M, = /Ms (A: ~p') and ;M„ = iM^ik ^ -k'). We have 

iM, = -/ „ ^ , {p-p'fG{po-p'o.P-p') (29) 



X [Av^pp' -P-P'] [Av^kk' - K ■ K'] , 



iMu = -i 



pp'-j 
8\/3/ 

,.2„;' D i-'l r^..2 n' 



(/9-^')'G(po-^(„p-fc') 



IV. VARIATIONAL CALCULATION 

We are now in a position to compute the viscosity due to 
binary collisions. We have to solve the linearized Boltzmann 
equation Eq. ( fTSl l with the scattering amplitude determined 
in the previous section, and then compute the viscosity us- 
ing either Eq. ^ or Eq. ( fTSI ). This task is simplified by a 
number of useful properties of the linearized collision oper- 
ator —Fij[g{p)]. The collision operator is a linear operator 
on the space of functions g{p). With a suitably defined inner 
product this operator is hermitian and negative semi-definite. 
As a consequence it is possible to compute transport proper- 
ties using eigenfunction and variational methods 1231. 

We elect to use the trial functions 



8ip)=p"Y.bMp), 

1=0 



(30) 



where « is a parameter that we choose for best conver- 
gence ifsoll . The orthogonal polynomials Bs{p) of order s are 
defined such that the coefficient of the highest power p'' is 1 
and that the orthogonality conditions Isill 



(0) 

+ff^)PijP"Br{p)B.{p) -A„5„, (31) 



d'p f. 



are satisfied. Starting from Bq = I we can recursively deter- 
mine all the Bs{p). This also defines the normalization factors 
A„. The polynomials Bs{p) are a generalization of the Sonine 
(modified Legendre) polynomials to Bose-Einstein statistics 
and linear dispersion relations. 



Inserting the trial function into Eq. ([7]) we find the following 

expression for the viscosity 



^Isip)] 



2v2 
5T 



.s=0 



fl^\l+fj,') (32) 



(271) 32£ 



X pijp"pijB,{p) 

2v2 ~ ^ 2v2 

-'^ .s=0 -'^ 

Alternatively, we can use the trial function in Eq. ( fTSI ). We get 

2 f d^p 



11 [^(P)] 



(2::) 



■ Pijg{p)piAg{p)] 



; £ b,bM.,, (33) 

s.t=0 



where Mj, are the matrix elements of the linearized collision 
operator 



(0) ^(0) 

k' 



X p"B,{p)pij [B,{p)fpij+B,{k)k"kij 

with the four-particle phase space factor 

^ _ d^p d^k d^k' d^p' 
P'"-^'''' ^ {2%)^2Ep {2%)^2Ek {2%)^2Ek> {2%)^2E,/ 

X {2Ky5^^\p + k~k'~p')\M\^. 
Eqs. (|32] | and (|33]| are consistent provided 



(34) 



(35) 



2v^ 

5r 



'sO- 



(36) 



This is a simple linear equation for /?, which is solved by 



bi 



2v^ 



-1 



' 





(37) 



Once bo is determined we can extract the viscosity from 
Eq. ( |32] |. In practice we pick a value for n and study conver- 
gence as the number of orthogonal polynomials is increased. 
What is nice about the method is that this is a variational pro- 
cedure. One can show that 11291] 



il> 



4y4 (fcoAoo)^ 
257-2 ^.jbsbtM,, 



(38) 



for any « and sets of b^. The condition that the bound is opti- 
mized with respect to the expansion coefficients bs is equiva- 
lent to the consistency condition Eq. ( |36] l. 

The scaling behavior of the viscosity with respect to the 
temperature and the Bertsch parameter t, is easily derived. We 
scale all momenta as p ^Tp/ v. Using Ep = v|p| this fixes the 
scaling of the energies. All terms in the scattering amplitude 
M have the same scaling behavior, except for a sub-leading 
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correction due to the self-energy insertion in the phonon prop- 
agator In terms of scaled momenta the phonon propagator 
G{pQ,p) can be written as 



iG{po,p) 



1 



T^pl-p^ + i^VHT/f^rimtipo^p) 



(39) 



The scaling of the scattering cross section is |M|^ ~ 
^^v^{T //u)^, and the self energy term induces corrections that 
are functions of ^^/^{T //j)"^ . We find 

Aoo =-jj^Aoo, (40) 

where we have dropped the corrections due to the phonon 
self-energy. At the leading order in the polynomial expansion 
8ip)^p"bo: 



A a2 



> 



Av" A, 



00 



00 



(41) 



25r2Moo 25v3r5^3^oo 25 Mqo 



where we used /j = t,Tf. An interesting dimensionless quan- 
tity to consider is the ratio of viscosity r| to the entropy density 
s for comparison with the conjectured bound discussed in the 
introduction, Eq. The phonon gas entropy is 



from which we obtain 



n ^ 72 ,,Al, (Tp 



> 

s 5571-^ ^ Moo V T 



(42) 



(43) 



In the calculation of M^i in Eq. ( |33] | the phase space inte- 
gral can be reduced to a 5-dimensional integral. Of the origi- 
nal 12-dimensional integration variables four integrations are 
removed using the energy-momentum conserving delta func- 
tion 5^^' {p + k + p' + k'). We choose to constrain the three- 
momentum p' and the magnitude |fc'|. Three more integra- 
tions can be removed as follows: without loss of generality we 
define the three-momentum p = pz along the z-axis eliminat- 
ing two angular integrations. Now, among the angular integra- 
tion variables only the z-axis projection of the three-momenta 
k (p - k = cos 9) and k' (p - k' = cos 9'), and the angular sepa- 
ration between k and k' [k- k' ^ cosy) are relevant. Thus the 
five remaining integration variables are: two magnitudes |p| 
and |fc|, two angles 9 and 9', and the angular difference (j) — (|)'. 
The 5-dimensional integration is done using the Monte Carlo 
routine Vegas ll32ll . 

In addition to varying the parameter n in the trial func- 
tion g{p) = p"T,sbsBs{p), we check for convergence as we 
increase the number of terms inside the summation. From nu- 
merical experiments with integer «, we find the maximal, con- 
vergent results for « = —1. In Fig.m we show {IQT //j)^r\/s 
at r = O.OOI/J with ^ = 0.4 for n — —2,-1. Convergence as 
the order of the polynomial used in the trial function is varied 
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FIG. 4: Numerical results for the scaled shear viscosity to entropy 
density ratio (10r/^)'^r|/i- at T = 0.001^ with ^ = 0.4 as a function 
of the polynomial order i'+ 1 of the polynomial We show 

results for two values of the variational parameter: square n = — 1, 
triangle n = —2. The straight lines connecting the numerical results 
are to guide the eyes. 
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FIG. 5: Shear viscosity to entropy density ratio as a function of 
T/Tf. The dots are numerical results, and the solid curve shows 
the power-law fit given in Eq. ( I44t . with t, = 0.4. 



is demonstrated. The n = —2 solution at leading order of the 
polynomial expansion starts small, and then converges to the 
« = — 1 result. This is expected since the « = — 2 trial function 
at second order of the polynomial Bi (p) ^ p contains the trial 
function with « = — 1 . 

Fig. |5] shows r[/s at three temperatures for the best trial 
function with n = —1. The data is very well described by 
the functional form 



-=7.7x 10 
s 



7^5 



(44) 



which is also shown in the figure. The numerical results in 
Fig.|5]are stable to about 1%. A comparison with the conjec- 
tured viscosity bound 1 / (47t) is shown in Fig.|6l The bound is 
violated for T > O.lTp, which is close to the measured critical 
temperature 7^ = (0.29 ± 0.02) ll28ll for superfluidity where 
the phonon calculation is not reliable. In the region where the 
phonon calculation is reliable, the viscosity bound is satisfied. 
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FIG. 6: Shear viscosity to entropy density ratio T|/j = 7.7 x 
10^^^^^ {Tf /T)^ compared to the proposed bound l/(47l). We show 
results for three values of the universal parameter ^: ^ = 0.4 (solid 
curve), ^ = 0.3,0.5 (short-dashed curves). The critical Tc = 0.29Tf 
is indicated. 



FIG. 7: Shear viscosity to entropy density ratio r\/s as a function of 
T/Tf. Solid curve: low temperature behavior of r\/s from Eq. J44t 
with ^ = 0.4, dashed curve: high temperature behavior of r|/i' from 
Eqs. l l45b and ( I46t . long-dashed curve: proposed viscosity bound 
l/(4n). Dots are data from jH. The critical Tc = 0.29rf is in- 
dicated. 



In Fig. |7]we compare our results to calculations in the high 
temperature limit and to experimental data. The high tem- 
perature results are taken from ifisll . These authors computed 
the viscosity due to binary fermion collisions. The free space 
cross section is proportional to 1 /k^. In the high temperature 
limit the infrared divergence is effectively cut off by the ther 
mal momentum (mT)^^^. For T 3> Tc 111311 



15 



3201 



(mr)3/2. 



In this limit the entropy density is that of a classical gas 



2V2 



3/2 



log 



4 tVz 



(45) 



(46) 



The data points are based on a hydrodynamic analysis 113 311 of 
experimental data on the damping of collective excitations in 
a unitarity Fermi gas |34t|. 

We observe that the naive extrapolation of the high T 3> 71 
and the low T ^ curves cross at around T w 0.27/7, which 
is indeed close to the transition temperature Tc ~ 0.297>. This 
crude extrapolation of the two limiting curves for r|/i suggests 
that the viscosity minimum is about a factor 5 above the vis- 
cosity bound. This is quite consistent with the experimental 
data. We also note that the experimental data show the ex- 
pected increase in r|/i for T ^ Tc, but they do not show the 
rise for T ^ 71 . This may be related to the fact that the phonon 
mean free path becomes so large that it is comparable to the 
size of the experimental Fermi gas sample and hydrodynamics 
does not apply. 



V. CONCLUSIONS 

We computed the shear viscosity of a cold unitarity gas in 
the superfluid phase. For T <C 71- ~ 7> the viscosity is domi- 
nated by phonons, and the leading order effective Lagrangian 



for the phonons is characterized by a single universal param- 
eter ^. This parameter can be extracted from the ground state 
energy of the unitarity gas. 

The calculation is based on the linearized Boltzmann equa- 
tion, and only the leading order 2^2 phonon scattering pro- 
cesses are included. The shear viscosity is determined using a 
variational procedure. We find that the shear viscosity scales 
as ri > 9.3 x \Q-%^T^ /{v^T^). This result can be combined 
with high temperature calculations of the shear viscosity to 
provide an estimate of the location and magnitude of the vis- 
cosity minimum. We find that the minimum value of r|/i oc- 
curs close to Tc, and that the value of r\/s is likely to exceed 
the proposed viscosity bound. A similar viscosity minimum 
is expected to occur in QCD. At low temperature the viscosity 
is dominated by weakly interacting Nambu-Goldstone bosons 
(pions and kaons) ||30[ [3ll [35ll . and at high temperature the 
viscosity is governed by weakly interacting quarks and glu- 
ons iH. 

There are a number of issues that deserve further study. The 
viscosity of the superfluid unitarity gas has the same 1 /T^ be- 
havior as the viscosity of liquid Helium at low temperature. 
In the case of liquid Helium the viscosity is believed to be 
dominated by 1 ^ 2 processes. On-shell phonon splitting pro- 
cesses can only occur if higher order corrections to the effec- 
tive Lagrangian lead to a concave phonon dispersion relation 
[Ep ~ v|p|(l +YP^) with Y> 0]. For the unitarity Fermi gas we 
do not know whether this is the case. It is known that the Bo- 
goliubov spectrum of a weakly interacting Rose gas has y > 0, 
so the role of these processes can be studied in an expansion 
around the Bose-Einstein limit. 

In general we would like to extend the calculation to higher 
temperatures. In the vicinity of the Tc we expect both bosonic 
and fermionic excitations to play a role. A possible start- 
ing point in this regime is provided by the expansion around 
d — A — z spatial dimensions proposed in [37!, 38]. It is also 
interesting to improve the high temperature calculations by in- 
cluding correlations between the fermions. Some steps in this 



direction were taken in 111311 . The result for time-like momenta and /?o > is 
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I-^l/'o^P) = d\k\ (A4) 



8v^fc^-8pov|fc|-3vV + 3pg)^(l+/£-/(:' , 



APPENDIX A: ANALYTIC FORM FOR THE SELF-ENERGY 

The imaginary part of E(po,p) arises from the pole terms 
in Eq. (|24|. We find (see [16]) 



ImE(/7o,p) 



3e3/2 



d k s\S2 



-H{P,K,si) 



16^/3/,,,^=±7 (27t)3 4£,£p_fc- 

b{po-s\Ek-S2Ep_k), (Al) 



1 , .(0) .(0) 



with 



H{P, K; s) = [pQ {2P-K-K^)+ ko{P^ -2P-K) 



ko=sEk 

(A2) 



There are four terms, corresponding to s\^S2 = ±1. Terms 
with s\ ^ S2 contribute for space-like momenta v|p| > pQ, 
and terms with s\ — S2 contribute for time-like momenta. For 
space-like momenta we get 



ImE(/7o,p) 



3V37t^V2p2 



128v|p|/ jEoz^ 



d\k\ 



(A3) 



(Sv^fc^ - 8pov|fc| - 3vV + 3pg)^ (C - /iiV , 



and ImX;(— /?o,p) = — ImE(/:'o,p). These integrals can be 
computed analytically in the limit of small momenta. In the 
space-like region we have \po\ < \'\p\ <C v|fe| ^ T. This im- 

Pl-/iS-/i|V.o-/'o(/,S)'and 



Sa* V \p\ 



For time-like momenta \po\ ^ v\k\ ^ T ^ v|p|, and 



lmnpo,p)^'-^^yVo 

exp(^) + l 
exp(^)-l 



exp(#) + l 
exp(#)-l 



0(-po) 



©(po) (A6) 
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